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Abstract 



The kinetic equations describing irreversible aggregation and the scahng approach de- 
veloped to describe them in the limit of large times and large sizes are tersely reviewed. 
^ , Next, a system is considered in which aggregates can only react with aggregates of 

' their own size. The existence of a scaling solution of the kinetic equations can then 

, be shown rigorously in the case in which the total mass of the system is conserved. A 

' large number of detailed properties of the solution, previously predicted by qualita- 

tive arguments, can be shown rigorously as well in this system. In the case in which 
\^ ' gelation occurs, some sketchy rigorous results are shown, and numerical evidence for 

, the existence of a scaling solution is presented. These are the first explicit examples 

' of typical scaling behaviour for systems exhibiting gelation. 

1 Introduction 

S ' In many systems one finds that irreversible aggregation of "clusters" A(m) of mass m oc- 

. curs and plays an important role. Particular instances are aerosol physics, vi^here suspended 

^ I particles coagulate due to van der Waals forces, polymer chemistry and astrophysics (where 

^ ' the clusters may be quite varied, going from galaxies to planetary systems). In all these 

systems one is among other things interested in the cluster size distribution as a function 
of time. To obtain it one then relies on kinetic equations, derived from the following 
assumptions: let the reaction 

A{m)+A{m) — > A{m + m') (1.1) 

K(m,m') 

occur at a rate K{m,m'), that is, let any two clusters of size m and m' react on a time 
scale given by K{m,m')^^ . If one further assumes that no spatial correlations build up, 
in other words, that we may use a mean- field model in which the probability of encounter 
of two clusters of masses mi and m2 is proportional to the product of the concentrations 
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of such clusters, one obtains the following kinetic equations for the concentrations c(m, t) 
of clusters of mass m: 



oo 



c(m,t) = - / dmi dm2K{'mi, 1112) €{7711,1)0(1712, t)x 







2 

X [6(771 — mi — m2) — 6(m — mi) — 5(m — m2)] ■ (1-2) 



These are an infinite set of coupled nonlinear ODE's, which represent a challenging prob- 
lem. Few exact solutions are known, which are reviewed in jl7j . 

From the form of the reaction scheme as well as from the equations (|1.2)) it is 

clear that at the formal level the total mass M 

M= 77ic(77l,t)d77l (1-3) 

Jo 

is conserved. It is well known, however, that this is not generally true. If K(7ni,7n2) 
grows fast enough with the mass, a finite amount of mass can escape to infinity in finite 
time. This phenomenon is known as gelation and is found explicitly for the reaction rates 
K(7ni,m2) = mims HI [H HI US] • 

Existence and uniqueness results are known for (|1.2j) . Let me state two which are 
representative: White [23] has shown that under the hypothesis 

K (rrii, 7712) < 0(7711 + 7712) (1-4) 

and if the initial condition 0(771, 0) has finite moments of arbitrary order, then the kinetic 
equations H1.2|) have a unique solution for all positive times and this solution preserves 
mass. A complementary result obtained in ^31 states the following: under the hypothesis 

K(mi,m2) < r(7ni)r(m2) 

r(m) = o(m), (1-5) 

the kinetic equations ()1.2|) have at least one (positive) solution for all positive times, which 
need not conserve mass. This was significantly extended in |1 H I12j. where it was shown 
that if 

K(mi,m2) = r(77ii)r(m2) + a(mi, 7712) (1-6) 
a(mi,m2) = a(m2,mi) < Ar(77ii)r(m2) (l-^) 

for some positive constant A, then a global positive solution exists for all times. In this 
case all requirements of growth on r(m) could be dropped, but the rate must now be 
dominated by a contribution of the product type. There also exists a rigorous result by 
Jeon |S], showing that no substantial improvement on White's theorem is possible: more 
precisely, for kernels satisfying (|1.5|) and additionally 

K(mi,m2) > e(mi77i2T (a > 1/2) (1.8) 

for some e > 0, it is shown that gelation always occurs. 

The rest of this paper is organized as follows: in Sections 2 and 3 I tersely review 
the definitions and predictions of scaling theory for non-gelling and gelling systems re- 
spectively. In Sectiond 4 and 5 I obtain specific results concerning the so-called diagonal 
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kernel in the scaling limit. In Ssection 4 I treat the non-gelling case, for which both the 
existence and many qualitative properties of the solution can be rigorously established. 
These results are summarized in Theorem^ In section 5 I discuss both the gelling case for 
which I am not able to show existence rigorously and the marginal case in which the ho- 
mogeneity degree is one which, although non-gelling by White's theorem, lies outside the 
domain of application of Theorem^ Existence and qualitative properties of the solution 
are summarized in Theorem |2j I present some conclusions in Section 6. 

2 Scaling Theory: The Non-geUing Case 

Beyond the above results on existence of solutions, there has always been a great interest 
in knowing the qualitative behaviour of the solutions either in the large-time limit or, for 
gelling systems, for times immediately before the singularity. This has been attempted by 
the so-called scaling description of the Smoluchowski equations, to which we now turn. 
Here we consider the non-gelling case. Scaling theory is then concerned with the large-size 
behaviour of the cluster size distribution c(m, t). It therefore deals exclusively with a limit 
in which time and size go jointly to infinity in a manner to be described shortly. 

To proceed, I first need to classify the rates K{mi,m2) according to the way in which 
they behave for large values of the arguments. One first assumes that they are at the very 
least asymptotically homogeneous, that is, there exists a A such that 



It turns out that only the homogeneous limit K{x, y) is relevant in the scaling limit. One 
then needs to classify homogeneous kernels K(x,y) further according to the way they 
behave when x and y are very different. One may write K{x,y) in the form x^k{y/x) 
and, since K{x,y) is symmetric, one has k{l/^) = The reaction kernel K{x,y) is 

therefore uniquely specified by the function on the unit interval, which may further 
be chosen arbitrarily. For simplicity, let us now assume that — > 1 as ^ — > 0, which 
defines the exponent ^. It then follows that k{S^)/^'^ — > 1 as ^ — > oo, where z/ = A — /x. 
The exponent v should satisfy u < 1 for the classical existence theorems to hold, for which 
H1.4|) is required. As stated above, there exist more general theorems 112| . extending 
to some cases for which u > 1, however I shall not consider these here. 

We now turn to the description of the scaling theory for non-gelling systems: it was 
originally developed by Friedlander [HJ [7j and later extended by Ernst and van Dongen 
|22j . The underlying idea is the following: if K(mi,m2) is exactly homogeneous of order 
A, then (|1.2j) may have self-similar solutions of the form 



where W is an unimportant normalization and the exponent z is given by z = 1/(1 — A). 
It is a straightforward computation to verify that this is the case if and only if the function 
$(x) satisfies the integral equation: 



K{x,y) = lim s K{sx,sy) 



(2.1) 



c{m,t) = Wt~'^^^{mf') 



(2.2) 




(2.3) 
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and the conditions 

x^{x)dx = 1 

^{x) > 0. (2.4) 

The problem of showing the existence of self-similar solutions is therefore reduced to the 
corresponding problem for (|2.Hj) . This, however, turns out to be quite difficult. Existence 
has recently been shown for three different specific (and rather typical) forms of the rate 
kernel K{x,y) in jSj, namely: 

Zi(x,y) = (x"+y")(x-^ + y-^) (a G [0, 1), /3 G [0, oo)) (2.5a) 

Z2(x,y) = (x"+y")^ (a,/3G [0,oo)) (2.5b) 

Z3(x,y) =xV + xV (a,/3G(0,l)). (2.5c) 

For these cases Fournier and Laurengot have shown existence (though not uniqueness) 
of a positive solution to 1)2. 3(1 decaying faster than every power for large argument and 
satisfying integrability conditions near the origin corresponding, for Ki and K-^, to the 
results previously conjectured in the literature and described below. The kernel K2 is in 
a sense, as we shall see shortly, marginal and results in this case are far more fragmentary. 
In very closely related work Escobedo et al. |^ show existence of scaling solutions for a 
kernel of the same type as if 3, but in the parameter range — I<a<0</9<1 and 
< A < 1. 

Let me describe what is supposed to hold for the behaviour of the solution <I>(x) of 
1)2. 3|) with the boundary condition 1)2. 4|) . For the behaviour near the origin it has been 
convincingly argued by power counting arguments (see for example |22j ) that the behaviour 
is different according to the sign of the exponent /i: if ^ > 0, then ^(x) behaves as x~^^^'^\ 
whereas if ^ < it behaves as a stretched exponential: 

const. ■ x~(^+^) (/i > 0) 

const. ■ exp const. ■ x"'^') (/i < 0), 

where in the second case a power-law behaviour presumably exists as a prefactor. Finally, 
if fi = 0, the behaviour near the origin is non-universal. This case is considerably more 
complex than the other two, and I shall not consider it any further. As an example, we 
may note that out of the kernels in (|2.5al2.5bl2.5cjl . Ki has jj, = —(3 < 0, K2 has = 
and K3 has ^ = max(a, (3) > 0. We therefore see that the three cases are represented 
in the results of 0. As noted there, it could be proved that, for Ki, <^{x) decays at the 
origin faster than any power, whereas for K-^ it could be shown that 

/ xP^{x)dx < 00 (2.7) 
Jo 

for all p > A = Q + /3, which is a precise form of the statement made in ()2.6|1 and therefore 
presumably optimal. 

The behaviour for large values of x has been found by similar arguments in The 
result is that <&(x) decays exponentially with a power-law correction, that is 

^{x) = const. ■ x~^ exp{— const. ■ x), (2-8) 

where this defines the exponent 9. One then finds 9 = X for i' < 1, whereas if = 1 the 
behaviour is non-universal and 9 cannot be evaluated in general. 



Hx) = \ (2-6) 
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3 Scaling Theory: The GeUing Case 

Let us now pass to the treatment of both the gelling case and that in which the degree of 
homogeneity A = 1. In the former, it is generally assumed, on the basis of numerical work 
and exactly solved models, that most of the mass is contained in aggregates of non-singular 
size. In formulae 

r-M 

lim lim / mc{m,t)dm = 1. (3-1) 

Since the distribution mc{m, t)dm is concentrated on small clusters and does not, therefore, 
display the singular behaviour taking place at t^, we should look at higher powers, in 
particular at the distribution m?c{m, t)dm. Let us assume pZ| that, after rescaling m by 
an appropriate factor, this distribution converges weakly to a distribution x'^(^{x)dx, that 
is, 



lim 

t— »oo 



/■oo 

Mi{t)-^ / dmm^c{m,t)f{m{tg-tf/'') 
Jo 



oo 







oo 



M2{t)= I m^c{m,t)dm. (3.2) 
Jo 

Here the rate at which the typical size diverges is given traditionally by the exponent — 1/cr. 
In order for these equations to be closed, however, we need further knowledge concerning 
the behaviour of M2{t) near tg. This can be obtained by the following argument. It is 
clear that for gelation to occur the second moment must diverge at tg. The concentrations 
c{m,tg) have therefore some kind of power-law behaviour which we may denote by m""^. 
If we now take M2 (t) to be approximated by 

/ m^c{m, t)dm = {tg - (3.3) 

Jo 

and if we additionally assume that the exponent r characterizing particles of fixed size 
m as t ^ tg coincides with a power-law singularity at the origin of the scaling regime, 
in other words, if we assume that particles of size very small with respect to the typical 
size have the same power-law singularity as particles of fixed size when that size goes to 
infinity, then we are led to the additional boundary condition 

= const. ■ x~'^ {x 0). (3.4) 

These assumptions from now on lead in a rather straightforward way to the following 
extension of (|2.3|) for gelling systems: 



dx [{3 - T)f{x) + xf'{x)] x'^^{x) = J dxdy K{x,y)x^{x)^{y)x 

x[{x + y)fix + y)-xf{x)] (3.5) 

for all continuously differentiable functions /(x). The problem is therefore reduced to a 
kind of non-linear eigenvalue problem: one must find a value of r such that there exists 
a solution $(x) of ()3.5() satisfying the boundary condition (|3.4() and decaying sufficiently 
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fast at infinity. We shall show in an explicit example, namely the so-called diagonal kernel 
K{x,y) given by x''^^^6{x — y), how the problem may be solved at least in practice and 
give strong numerical arguments for the existence and uniqueness of a solution. 

For the exponent a characterizing the divergence of the typical size, one obtains ^7] 

o- = 1 + A - r, (3.6) 

in stark contrast with the non-gelling case, in which the corresponding exponent z could 
be expressed directly in terms of the known value A. Here, on the other hand, 1)3. 6|1 
only yields a relation between two equally unknown quantities. There does exist jl7| an 
interesting relation between a and the correction to scaling exponent A of <I>(x) defined 
by 

^>(x) = const. ■ x-^ [1 + 0(x^)] . (3.7) 

It is then conjectured that (see [T7j for details) that 

A = CT = 1 + A-T. (3.8) 

In this paper we shall show rigorously for the case of the diagonal kernel that, if a solu- 
tion to the non-linear eigenvalue problem described above actually exists, then the above 
relation holds. 



4 The Diagonal Kernel: The Non-gelling Case 

In the following we study in greater detail the so-called diagonal kernel given by K{x, y) = 
x^^^5{x — y) Physically, this means that only aggregates of the same size react with each 
other. This might at first be viewed with suspicion: for example, if one starts with 
initial conditions of the type c(m, 0) = 5{m — mg) then a very singular solution arises, in 
which only sizes of the form 2^mQ occur. However, we may start from less singular initial 
conditions and should also remember that in the scaling limit any kernel which only allows 
reactions between nearby sizes would be described asymptotically by this kernel. 

It follows from White's theorem that A < 1 is non-gelling, whereas a specific study 
of this particular model shows that the system undergoes gelation for A > 1^. In 

this section we only consider A < 1, as the case A = 1 turns out to be highly singular and 
best treated by methods akin to those used for the gelling case. One can then be quite 
explicit about the properties of the solution of (|2.3j) : 

Theorem 1. Let the reaction kernel K{x,y) be given by x'^^^(5(x — y) where A < 1 and 
let A be the unique positive solution of the transcendental equation 

1 + A i-2(^-mA+i) 

Then there is an e > and a solution <I>(x) of \2. such that 

1-A 



$(x) = X 



-(l+A) 



1 - 2^-1 



cx 



+ O (x^+^) 



(4.2) 



^Note that this does not follow directly from Jeon's more general result in 
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Here c is an arbitrary positive constant, which can be varied so as to set the normaliza- 
tion \2.4\j - Furthermore, if^{x) is of the form \4-'^ , then it is uniquely determined by c. 
Finally, x^^^'^{x) is monotonically decreasing on [0, oo) as well as bounded by an expo- 
nential function C exp(— ax) for some positive constants C and a. It is also C°° on [0, oo) 
in the variable x^. 

Remark. Note that the Theorem makes no claim concerning uniqueness of ^(x) if <I>(x) 
does not have the form (|4.2|) . In particular, solutions having divergent second moments, 
as found in the case of the additive kernel by Menon and Pego |2J are not excluded for 
this kernel. It would certainly be of great interest to know whether such solutions exist. 
Finally, x^^'^^{x) can be shown to be analytic in x^ in a neighbourhood of the origin. I 
believe, but cannot prove, that this is actually true for the entire positive real axis. 

Proof. Let us first write ()2.3() out explicitly in the case of the diagonal kernel. One finds 



= / dyy^+^<^{yf 

Jx/2 



t:/2 

which can easily be rewritten as a "non-local" differential equation: 



(4.3) 



A = ^A+2^(^)2 - i (I) (I) ^ (4.4) 

We now perform the following change of variables 

y = —y (4.5) 

which casts (j4.4j) in the form 

(1 - A)-^ [y^Piy)] = - 2^-V(2^^iy)^ (4.6) 

dy 

It is clear that the only way in which such a function can tend to a constant value tpQ at 
the origin arises if 

^0 = 3^^. (4.7) 
We therefore begin by searching for solutions of the type 

^(y) =Vo + '/>(y). (4.8) 
(j){y) then satisfies the equation 

(1 - A) [ym]' = 2V^o [m - V(2'-'y)] + <i>{yf - 2^'^ <t>{2^-^ y)\ (4.9) 

which always has the trivial solution (j){y) = corresponding to the solution of 1)4. 3(1 given 
by $(x) = iI!qx~^'^^^\ The existence of such a formal solution, which cannot satisfy the 
normalization condition (|2.4I) is a well-known fact for all kernels of type I. We should 
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therefore look for ways of violating uniqueness for the solution of (|4.9)) , which can clearly 
only happen near the origin. 

If we wish to construct a solution starting at and having the form (|4.8)) . then 4>{y) 
must initially be small. This suggests looking at the linearization of ()4.9|) . which can be 
solved by an appropriate power law y^. One then obtains the transcendental equation 

Here the second equality defines F{A). The solution is found to exist and to be unique: 
indeed, one finds F(0) = 1 > 1/2, but F{A) becomes less than (1 + A)/2 as A — > ±00, so 
1)4.10(1 must have both a positive and a negative root. But since F{A) is concave, there 
are no more than two roots, so that the positive solution is unique as stated. Further, 
since F(A) > 1 for all A > 0, it follows from (limH) that A > 1. 

It now remains to show that a solution (p{y) of the type described above indeed exists 
in a sufficiently small neighbourhood of the origin. To this end we introduce 

<A(y) = y^[-c + %)], (4.11) 

where c is an arbitrary constant, which we shall later take to be positive, but its sign is of 
no importance for the local result I am now proving. Using the transcendental equation 
(BTTnl) satisfied by A, we may now recast the equation ()4.3|) as follows 

= (T^TyTA dww'^{2i^oh{w)+w''[-c + h{w)f]. (4.12) 

To show that this fixed point problem has a unique solution, it is sufficient to find a closed 
subset of a Banach space which is left invariant by the operator 

nhKy)= / _ dww''{2^Pohiw)+w^[-c + hiw)f] (4.13) 

and in which T is a contraction. Define Coo,e(y) to be the space of continuous functions 
vanishing at the origin and so that there exists a C such that 

C-^5(C) <C {0<C<y). (4.14) 

If one then defines 

\\h\\oo,e= sup \w-^h{w)\, (4.15) 

0<w<y 

it is easy to see that the operator T maps Coo,e{y) on itself as long as e < A. It is further 
straightforward to show that 

||T[/i]|L,, < C{e) \\h\\^^^ + a{y) + %), (4.16) 

where C(e) < 1 whereas a{y) and b{y) go to zero as y ^ 0. From this follows that T leaves 
a ball of appropriately small radius R invariant if e is fixed and y is small enough. 
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To show that T is in fact a contraction on such a ball for y sufficiently small, one 
computes the norm of 

T[hi]-T[h2] = - 1^^-^ r dww^{2xPo[hi{w)-h2(w)] + 

[hi{w) - h2iw)] [hiiw) + h2{w) - 2c] }. (4.17) 
For the first term one shows 

^ dww^2ij}Q[hi{w) — h2{w)] < — [ dww^^'^ 'Whi — h2\\^ ^ ■ 



A-1 



(4.18) 



It then readily follows from the equation (|4.1fl)) defining A that the numerical constant 
before ^ is strictly less than one if e > 0. For the second term, on the other hand, 

one finds, as is easily checked, that it is bounded by 

c^(y)(ll^illoo,e + ll/^2lL,J 11/^1 -/^2lL,. (4.19) 

where C{y) goes to zero as y does. Since we have restricted ourseleves to a ball of given 
radius R, the prefactor ^ + ||^2|loo e bounded by 2i?, so that the entire mapping 

T is indeed a contraction for y small enough. 

To show global existence and decay to zero, we must first choose c > 0. From this 
immediately follows that there is a neighbourhood of zero in which the solution V(y) is 
monotonically decreasing. One now shows using ()4.6() that this property never disappears 
for any positive value of y. Indeed, if ^{y) were to increase for some value of y, it would 
have had to pass through a minimum at some value of y. Let yo be the first non-zero value 
of y for which ip^y) = 0. One then finds 

= V(yo)' - (1 - A)V(yo) - 2^-V(2^~'yo)' 

< (1 - 2^-')i^{yof - (1 - A)V(yo) (4.20) 



< ij{yo) (1 - 2^-^)Vo - (1 - A) 



0, 



which is absurd. Here the first inequality follows from the monotonicity of ip{y) before yo, 
and so does the second, since V'o is the initial value of Tp{y). Both are strict since ip'{y) ^ 
for y < yo. 

Once the local solution near the origin has been found, however, the equation H4.6|l 
can be considered as an ODE and the standard existence and uniqueness results apply: 
indeed, the function ■ip{2^~^y) is already known at y. In other words, if ipn{y) is defined to 
be '4>{y) restricted to the interval In given by [2(-'^~''^)("'~-^) , 2'^^"''^-'"'], then equation (|4.6|) can 
be viewed as an ODE connecting tp'^ and ipn-i- This then allows a recursive application of 
the usual existence and uniqueness results, since there is, by the local result just proved, a 
unique solution for n sufficiently negative. Furthermore, the integral form of the equation 
for ip{y), which is given by 

(1-A)y^(2/) = r dwiPiwf (4.21) 



2A-1 



y 
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guarantees the positivity of 4^{y) wherever it exists. With these remarks, standard theo- 
rems on ordinary differential equations yield global existence of ip{y) on the positive real 
axis. 

Proving that V(y) decays to zero exponentially fast as y ^ oo is now straightforward: 
indeed, from (|4.21|) and the monotonicity properties of 4'{y) follows immediately 

^(2/) < ^^^^('^^''y)' = V'o V(2'-^y)^ (4.22) 

from which one obtains exponential decay in the following way: since ip{y) < ipo for y > 0, 
by applying this identity repeatedly one eventually shows that there is an interval /„(, such 
that ip{y) < M < 1 on it. Now, if a; G /„, it follows that the inequality 1)4. 22p can be 
applied iteratively n — uq times to yield 

^{y) < V'o"^""""V(2("-"«)(^~')y)'"""° < V'o"^""""^ > (4-23) 

from which exponential decay in y follows, since y £ In means that (1 — A)(n — 1) In 2 < 
Iny < (1 — A)nln2. Note that the prefactor increasing exponentially in n corresponds 
only to a power law in y and is therefore a subdominant correction. 

I mention here an easier way of proving a local existence result, but which may be 
difficult to extend to an equally strong uniqueness result: one can make the Ansatz 

oo 

V^(y) =Vo + 5]c„y'^^ (4.24) 

n=l 

and plug this in (|4.(i|) keeping at first A a free parameter. A recursion for the c„ is easily 
derived, which is, however, in general inconsistent for ci. The consistency condition is, 
of course, the transcendental equation (|4.1()j) for A. ci can then be chosen arbitrarily, 
thereby determining the Cn uniquely. Convergence is also readily verified by inspection of 
the recurrence. This result shows as an added bonus that ^(y) is analytic in y^ around 
the origin. 

From the result shown we have therefore proved that for each negative value of ci there 
exists a solution of the original equation ()4.3|) of the form 



(4.25) 



for small values of x and that this function exists on the entire positive real 
positive function which decays exponentially to zero. It therefore satisfies the scaling 
equation and all the boundary conditions imposed upon it, since it is easy to see that the 
normalization to one of the first moment of $(x) can be obtained through an appropriate 
choice of ci , as the first moment is finite for functions of the type described by (|4.25|) . ■ 

5 The Diagonal Kernel: The Gelling Case 

Let us now consider the same diagonal kernel x^~^^6{x — y) but now for the case A > 1. 
We first consider A > 1 and then say a few words concerning the case for which A = 1, 
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which does not display gelation but for which the ordinary scaling ansatz described in the 
previous section is nevertheless inapplicable. 

The relevant equation (|H.5j) reduces for the diagonal case to 

4- Ix^Hx)] - (3 - t)x^^(x) = - 2-(^+3)x^+3$(x/2)2, (5.1) 

ax '- 

where ^>(x) should satisfy the boundary condition at zero 

^>(x) = const. ■ x~^ {x 0) (5.2) 
as well as a condition of rapid decay as x — > oo. One now introduces new variables 

^ = x^$(x) 

l + X — T 

C=T^T . (5.3) 

l + X — T 

Note that in our new variables the condition 1)5. 2|) translates into the requirement that ip 
be regular at the origin. We have also made the implicit assumption that 1 + A — r > 0: 
this is quite unproblematic, since through the relation (|3.6j) determining a it follows that 
this must always hold if we wish to describe a gelling system. 
The relation (|5.1() now transforms to 

^ = m'-2'-'-'i^{2-'-'C)'. (5.4) 
Local existence is here easy to show: one makes an ansatz in power series 

oo 

m=^akC''- (5.5) 
fe=o 

It is then straightforward to determine a recursion for the and to prove that it converges 
in a finite interval. No inconsistencies arise in this recursion, so we may simply set ag = 1. 
It follows that 

ai = 1 - 22^-^-3_ (5_g^ 

Let us then always limit ourselves to the case r > (A + 3)/2, so as to have an initially 
monotonically decaying function. Note here that this shows the claim made in the text, 
that the correction to scaling exponent of $(x) is 1 + A — r: indeed, we have shown that 
around C = 0) V'(C) is analytic, which together with the definition of C in (|5.3j) shows the 
result. In fact, the result obtained is stronger: it states that $(x) is of the form 

$(x) = x-^f , (5.7) 

where f{C) is a function analytic around zero. 

Since r > (A + 3)/2 one can show that this leads, as in the preceding section, to 
global existence unless the solution becomes negative. It is a crucial difference with the 
previous case that this can actually happen for given values of r, as we shall see shortly. 
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Figure 1. A solution of (|5.4I) corresponding to A = 2 and t = 2.6 plotted as a function of (: the 
dotted curve corresponds to 1/^ behaviour. Note the overall parallelism as well as the logarithmic 
oscillations of the function. The occurrence of such oscillations is intriguing, as they have been 
tentatively reported by Lee JHl for the scaling function in non-gelling systems of type 1. The curve 
shown here, on the other hand, is an inadmissible solution for a gelling system. 



One proceeds as before: let Co be the first value of C for which ip'{C) = 0. Since V'(C) is 
monotonically decreasing for C < Coj one finds using (|5.4|) that ip'{Co) < 0, which is clearly 
absurd. 

Performing a detailed qualitative analysis of (|5.4|) is much more difficult than the cor- 
responding task in the non-gelling case. I shall therefore appeal to numerical evidence, 
which seems to me quite compelling. Let us disregard considerations of uniqueness and 
look for a solution in the region r > (A-|-3)/2, for which I know that the solution is always 
decreasing. ip{C) therefore either crosses to negative values at some value of in which 
case it is surely not a physically relevant solution, or else it exists as a positive function 
for all positive values of C. In this latter case it is also easy to check that it must go to 
zero, as it cannot saturate to a constant non-zero value unless r = (A-|-3)/2. If r is barely 
larger than (A -|- 3)/2, it is relatively easy to see that ip{C) must cross to negative values 
and numerical work amply confirms this expectation. This is found by comparing (|5.4|1 
with 



f = /(C)'-(i + 6)/(o)2, 



(5.8) 



which is close, as for such values of r V'(C) is nearly constant. It then suffices to show that, 
for small values of e, the solution of (|5.8|) crosses the negative axis for values of C at which 
the value '(/'(2^^^^^C) is still close to its original value. Another asymptotic behaviour 
which is consistent with ()5.4() is a decay. Such a decay can in fact be shown to occur 
in the case r = A -|- 1, since ()5.4() can then be solved exactly and is found to have this 
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Figure 2. Two solutions of H5.4|l for A = 2 corresponding to two nearby values of r plotted as 
a function of x, that is ('1/(1+^^'^): for the lower curve r = 2.5542887 whereas for the upper one 
r = 2.5542888. Note how one crosses over to negative values, whereas the other seems to saturate. 
Note also the rapid decay common to both curves, which suggests the existence of an exponentially 
decaying solution of (|5.4() for some critical r. 

behaviour. Again, as shown in Fig.^for the case A = 2 and r = 2.6, such a decay is indeed 
observed for values of r somewhat larger than those for which a crossing to negative values 
was observed. What we want, however, is neither: negative values are unphysical and the 
is too slow to be acceptable. One needs a function having a rapid decay for large C- 
We may therefore ask whether between the two regimes there is a (conceivably unique) 
value of T for which this occurs. In fact, we anticipate exponential decay in x and hence a 
corresponding streched exponential decay in Fig. |21 gives a striking confirmation of this 
suggestion: it plots for A = 2 two quite nearby values of r. For the lesser, we see a crossing 
to negative values, for the higher, a saturation apparently to a constant (but theory states 
this to be impossible, so it is presumably a behaviour), but both exhibit over a very 
large range exactly the expected exponential decay: for clarity I have plotted V'(C) ^ ^ 
function of x on a semilogarithmic scale, for which the expected result is a straight line. 

Let us end by some remarks concerning the limiting case A = 1. If we refer to the 
results of the previous section for the non-gelling case, we find that they do not apply, as 
all equations become highly singular and the proofs all fail. The reason for this is to be 
found in the fact that if Theorem ^ were valid, $(x) would have an singularity at the 
origin, thereby precluding the use of the scaling hypothesis (|2.3|) . Indeed, for (|2.3|) to hold, 
the normalization condition (|2.4|) is necessary, and this is contradicted by the existence of 
such a singularity. We therefore treat the subject as if it were a gelling case and study the 
second moment as in (|3.2() . This leads ^7] to a typical size s(i) of the form 



s{t) = const. ■ exp{K\/t) 



(5.9) 
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and a function <I>(x) satisfying 
d 



dx 



x^^{x) = x''<i>{xy - {x/2y<^{x/2y 



(5.10) 



which is merely ()5.1() in the special case in which A = 1 and r = 2. 
We now state our results for this case: 

Theorem 2. For each oi < i5.1(J\) has a unique solution of the form 



= x 



1 + y^gfc^ 



k=l 



(5.11) 



x'^^{x) is monotonically decreasing on the positive real axis and hounded by an exponential 
function C exp{—ax) for some positive constants C and a. It is also Coo in [0,oo). 

Remark. Note that there are other solutions to (|5.1U|) which are not of the form (|5.11jl . 
In particular, the same proof as below shows that there are solutions of the form 



^0 + ^ akx 



kA 



k=l 



for all ipQ > where A is the unique positive solution of 
A 



2(1-2-^) 



■00 



(5.12) 



(5.13) 



It is readily seen, however, that these solutions bring no new scaling solutions, since the 
exponent A can be absorbed in the constant K of ()5.9|) . I also am not sure whether 
solutions with power-law tails having a form different from (|5.12|) could exist. 

Proof. The change of independent variable in (|5.3|) now becomes singular, so we merely 
change to 



tp = x^<I>(x). 



(5.14) 



Note that we still have the boundary condition ^{x) regular at the origin, corresponding 
to (|^ for T = 2. This then yields 



i^ixjlf. 

One finds a solution as follows: the ansatz 



(5.15) 



V'(x) = 1 + ^ a^x^ 



(5.16) 



k=\ 



leads, as is readily verified, to a recurrence for the which is consistent throughout 
and for which a\ can be chosen arbitrarily, after which all are uniquely determined. 
Convergence of the series (|5.16p in a finite interval around zero is also easily checked. Again 
choose a\ to be negative. It follows immediately that 'i/'(x) is monotonically decreasing in a 
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neighbourhood of the origin and hence, by an argument entirely similar to that presented 
in Section 4, wherever it exists on the positive real axis. The equation can also be rewritten 
in integral form as 



Jx/2 y 

which shows that ip{x) is always positive on the positive real axis. This leads, using the 
same approaches as previously, to a global existence result for the solution of 1)5. 10() . From 
(|5.17|) one also immediately derives the inequality 



from which, as in the previous section, exponential decay of 'tp{x) is readily deduced. ■ 

6 Conclusions and Outlook 

Summarizing, we have shown that the equations describing the scaling function for the 
Smoluchowski equations describing irreversible aggregation have a unique solution within 
a given family of physically reasonable solutions when the reaction rates are of such a 
nature as only to allow reactions between aggregates of identical sizes. Since this kernel 
is dominated by diagonal reactions, it belongs to the so-called type I kernels. For these 
the very existence of a solution to the scaling equation (|2.3|) has sometimes been doubted, 
both on theoretical and numerical grounds: theoretically, because of the existence of a 
non-normalizable solution of the type x~^^^'^\ which indicates that the constraint of fast 
decay at infinity must play an essential role in determining the solution. Numerical work 
on some kernels of type I also led to doubts concerning scaling jTUj, though these were 
later contradicted by who showed that the previous finding may have been due to 
non-monotonic behaviour of the scaling function for the rates under consideration. For the 
diagonal case under study, I have not only rigorously shown existence in the non-gelling 
case, but also monotonicity. For the gelling case, on the other hand, it was suggested by 
numerical evidence that whenever t > {X + 3)/2, there are only two possible behaviours 
in general, namely becoming negative or decaying algebraically to zero from above. Just 
at the border between these two unphysical behaviours, numerical evidence again strongly 
suggests that a positive and exponentially decaying solution exists. This then allows a 
determination of the exponent r. As stated already in ^7], it is not to be expected that 
the "standard" value r = (A -|- 3)/2 broadly quoted in the literature should be generally 
valid. In fact, it follows from the considerations made in this paper that no solutions 
of this type can exist for the diagonal kernel. For the limiting case A = 1, we have 
shown rigorously the existence of a solution of the scaling equation ()3.5() derived from 
an assumed weak convergence of the (appropriately rescaled) measure m?c{m, t)dm to a 
limiting distribution x'^^{x)dx. These are related to the modified scaling ansatz introduced 
by van Dongen |24j . No such explicit solutions have, to my knowledge, been derived before. 

Finally, it must be remembered that none of all this has any relevance unless it is 
somehow shown that the solutions of the time-dependent problem corresponding to the 
original Smoluchowski equations (|1.2|) can be shown to approach the scaling limit. It is 
possible that the recursive structure underlying the diagonal kernel might make this result 
somewhat easier to show. It remains, however, a formidable task, left to future research. 




(5.17) 




(5.18) 
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